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ABSTRACT 

Recent observations of UV-/optically selected, massive star forming galaxies at z ~ 2 indicate that the bary- 
onic mass assembly and star formation history is dominated by continuous rapid accretion of gas and internal 
secular evolution, rather than by major mergers. We use the Millennium Simulation to build new halo merger 
trees, and extract halo merger fractions and mass accretion rates. We find that even for halos not undergoing 
major mergers the mass accretion rates are plausibly sufficient to account for the high star formation rates 
observed in z ~ 2 disks. On the other hand, the fraction of major mergers in the Millennium Simulation is suffi- 
cient to account for the number counts of submillimeter galaxies (SMGs), in support of observational evidence 
that these are major mergers. When following the fate of these two populations in the Millennium Simulation to 
z = 0, we find that subsequent mergers are not frequent enough to convert all z ~ 2 turbulent disks into elliptical 
galaxies at z = 0. Similarly, mergers cannot transform the compact SMGs/red sequence galaxies at z ~ 2 into 
observed massive cluster ellipticals at z = 0. We argue therefore, that secular and internal evolution must play 
an important role in the evolution of a significant fraction of z ~ 2 UV-/optically and submillimeter selected 
galaxy populations. 

Subject headings: galaxies: formation — galaxies: evolution — galaxies: high-redshift — cosmology: dark 
matter 



1. INTRODUCTION 

In the cold dar k matter model of hiera rc hical struc- 
ture formation dBlumenfhal et all Il984t iDavis et all 
1985; SpringeletalJl2006) mergers are believed to play 
an important role in gal axy formation and evolution 
( Steinmetz & Navarro! 120 02). Mergers induce starbursts 
(Hernquist & Mihos 1995) and transform galactic morphol- 
ogy dNaab & Burkertll2003l) . Major mergers may drive the 
buildu p of the red sequence (Toomre 119771: iHopkins et al.1 
2008a). Dark matter models and many observ ations show that 
mergers are more freq uent at high redshift ( Fak houri & Mai 
200alConselicdl2003|) . 

However, there is growing evidence that a smoother growth 
mode may be important for the baryonic mass assembly and 
star formation history at high redshift. For example, the 
tight correlation between star formation rate (SFR) and stel- 
lar mass in UV-/optically selected star forming galaxies is 
indicative of buildup by c ontinuous gas inflow (Dad di et all 
I2007t iNoeske et al.1 l20Q7h. As part of the SINS survey 
dForster Schreiber et aD |2006; N. M. Forster Schreiber et 
al. 2008 in preparation), integral field spectroscopy of more 
than 50 UWoptically selected z ~ 2 star forming galaxies 
show a preponderance of thick gas-rich rotating disks and 
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only a minority of majo r merg e rs dForster Schreiber et all 
120061 iGenzel et alJ 120061 120081 IShaoiro et al] 120081) . In 
contrast, SMGs are probably short-lived maxim um- starburst 
galaxies undergoing dissipative major mergers dTacconi et all 
l2006tlBouche et alj|2007tlTacconi et"aTJl2008l) . Table Q] sum- 
marises key properties of these z ~ 2 galaxy samples. 

How do these observations fit into the concordance cosmo- 
logical model? Modern simulations of dark matter structure 
formation are robust and fixed by the cosmological param- 
eters. However, complicated baryonic physics makes it dif- 
ficult to model the evolution of galaxies and rep roduce, for 
exam ple, the high SFRs of these z ~ 2 galaxies (Daddi et al. 
120071) . 

Galaxies at z ~ 2 differ significantly from local galaxies. 
The central mass densities o f SMGs and of massive qu iescent 
galaxies at the same redshift dvan Dokkum et alj200 8 and ref- 
erences therein) are an order of magnitude greate r than those 
of local spheroids and disks (Tacconi et al. 2008). Also, the 
Z « 2 rotating disks are thick and turbulent, unlike local disk 
galaxies. These differences raise the question: what are the 
local Universe descendants of these z~2 galaxies? 

In this paper w e use the cosmolog ical dark matter Millen- 
nium Simulation ( Spring el et al.ll2005l $2) to investigate the 
possible role of major mergers in galaxy formation at z ~ 2 
(§3), and to consider the evolution of the z ~ 2 galaxies to 
z = 0(®. 

2. analysis of the millennium simulation 

2.1. The Millennium Simulation and its merger trees 

The Millennium Simulation is a cosmological N-body sim- 
ulation. It follows 2160 3 dark matter particles of mass 8.6 x 
10 8 /j _1 M Q in a box of 500/T 1 Mpc on a side from z = 127 to 
z = 0. There are 64 output times ("snapshots"), at « 300 Myr 
intervals at z < 3. The cosmology is ACDM, with fi m = 0.25, 
Q A = 0.75, Q b = 0.045, h = 0.73, n = 1 and cr 8 = 0.9. 
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TABLE 1 

Properties of galaxy samples at z « 2 



Galaxy 


SFR 


Halo 


Comoving number 


Major 


sample 


[Meyr 1 ] 


mass 


density 


merger 






[M ] 


[fc^MpcT 3 ] 


fraction 


SINS 


M 30-300 (ll) 


, n 11.84 3 x 

iu v 200 x 


1-2.2 X l(r 4(c - d) 


« 0.3 <e) 


SMGs (0 


k»750±300 ( s) 


jl+z r l.5 h -l (b) 
1 3.2 - 1 0.7 


1-2 X 1(T 5(C > 


w 1 <C) 



NOTE. — (a) N. M. For ster Schreiber et al. 200 8 in preparation ; (b) 
IForster Schreiber et al. 2006 ; (c) Tacconi et al. 2008 a nd references th erein 
; (d) BX/BM & sBzK galaxies with K < 20 : (e) IShapiro"eraT]|2008l ; (f) 
S(850/zm) > 5mjy ; (g) Genzel, priv. comm. 



In the Millennium Simulation, structures are identified in 
two steps. First , the Friends-of-Friends (FOF) algorithm 
dDavis et alJI 19851) creates a catalogue of FOF groups at each 
snapshot. The FOF groups represent dark matter halos. Sec- 
ond, bound s ubstructures are iden tified inside the FOF groups 
(SUBFIND; ISpringel et al.ll2001l) . so that each halo contains 
at least one subhalo. The Millennium merger trees are con- 
structed by finding a single descendant for each subhalo in the 
following snapshot, while the FOF groups themselves play no 
role in constructing the merger trees. 

In traditional merger trees, mergers are instantaneous, i.e. 
there is no information on their durations. Therefore, the Mil- 
lennium public merger trees 6 can be used to determine the 
merger rate, whi ch is merely a count o f the number of merg- 
ers per unit time dFakhouri & Mal2008l S. Genel et al. 2008 in 
preparation). However, they cannot be used to determine im- 
portant quantities such as the major merger fraction, defined 
as the fraction of halos undergoing major mergers at a given 
time, or the mass growth rate of each halo. 

2.2. Constructing new merger trees 

To derive merger fractions and mass growth rates, we must 
consider the finite physical durations of mergers. Therefore, 
start and end points must be defined. Also, to derive these 
quantities for entire dark matter halos rather than for subha- 
los, new trees have to be constructed based on FOF groups. In 
our procedure, the main subhalo in each FOF group is iden- 
tified and is then followed to its subhalo descendant, using 
the original subhalo-based trees. The FOF group to which the 
subhalo descendant belongs is defined as the FOF group de- 
scendant of the original FOF group. Thus, in our new trees 
each node is an entire FOF group, rather than a subhalo. If 
two subhalos merge while within a single FOF group we do 
not count this as a merger event. 

We identify a merger whenever two or more FOF groups 
at snapshot s have a common descendant at snapshot s + 1 . 
However, at this time the halos are not necessarily already 
physically merging, since they may still be well separated. To 
account for that, we track the distances between the subhalo 
descendants of the main subhalos of the original FOF groups. 
These subhalos represent, approximately, the centers of the 
entire halos. We then define the start point of the merger as the 
last snapshot at which this distance is still larger than the sum 
of the virial radii of the original halos (FOF groups). In some 
mergers the subhalos disappear before the distances become 
smaller than the sum of the original virial radii. When this 
happens the start point is defined as just one snapshot prior to 

6 The Structure catalogues and derived merger trees have been made public 
by the Virgo Consortium: http://www.mpa-garching.mpg.de/millennium 



the point where the subhalos merge and/or disappear. 

After a merger begins, one of the halos becomes a substruc- 
ture within the other. This substructure typically dissolves too 
quickly to be followed until the merger is physically complete. 
To overcome this problem we first estimate the duration of 
mergers r mergel , and define their end point as larger after the 
start point. To e stimate the durations, we co nsid ered the fit- 
ting fu nctions of lBovlan-Kolchin et al.l d2008) and lJiang et al.l 
(2008), which are based on simulations of mergers. Our 
results are qualitatively robust with respect to this choice. 
We present quantitative resu lts based on the orbit-averaged 
iBoylan-Kolchin et al.l (120081) fitting function for the dynami- 
cal friction merger time: T melger = 0.05 fa ' (1+r) g^, where r is 
the mass ratio and H(z) is the Hubble constant at redshift z- 

The accretion rate we associate with each merger equals 
the amount of accreted mass divided by the merger duration 
T 

1 merger- 

To summarise, we construct new FOF group-based merger 
trees, where each FOF group also holds information about in- 
ternal on-going mergers. The accretion rates associated with 
those mergers are summed up to obtain the total accretion rate 
onto the FOF group in question. The merger mass ratio is de- 
termined by the masses of the FOF groups at snapshot s, just 
prior to the appearance of a common FOF group descendant 
at snapshot 5+1. We define major mergers as those with mass 
ratios between 3 : 1 and 1 : 1 (with 1 : 1 being the most "in- 
tense" type of merger). If the most intense merger associated 
with a halo lies between 3 : 1 and 1:1, the halo is labelled as 
undergoing a major merger. 

3. GALAXIES AT z fa 2 

Fig. Q] shows halo number densities (shaded contours) and 
major merger fractions (red contours) as functions of halo 
mass and dark matter accretion rate for z ~ 2.2. It shows that 
the major merger fraction is an increasing function of specific 

dark matter accretion rate (ff^-), as both quantities increase 
towards the upper-left direction of the plane. This trend holds 
at all redshifts. 

The mean accretion rate scales with halo mass and redshift 

as 

(M DM ) » 35 M yf'Cl +zf 2 M\f (1) 

(where M\i = lQl f M )■ The lcr scatter equals w (Mdm) x 

(ffl) 02 > which reflects more the upwards scatter, although 
negative accretion rates do exist for some halos (because 
of tidal stripping or fluctuations related to the FOF algo- 
rithm). Our numerical results are in good agreement with 
the analytic approxim ation for the accretion rate presented by 
iNeistein et al.l 0006), wh ich is based on the extended Press- 
Schechter (EPS) m odel dPress & Schechtedll974tlBond et al l 
ll99lUBowerll991l) . The llSfeistein et all (120061) approximation 
has a somewhat stronger mass and redshift dependence com- 
pared to our results. For example, at z = and M = 10 12 M©, 
their accretion rate is w 10% higher, and at z ~ 3 and M = 
10 14 Mq it is a factor of w 2 higher. For halos of partic- 
ula r interest for this pa per, i.e. of M sa 10 12 M Q at z ~ 2, 
the Neistei net al.l (|2006) approximation exceeds our eq. [T|by 
w 30%. 

To compare our results to observed galaxies, we assume 
that the galaxies are the central galaxies of their halos. We 
convert dark matter accretion rate (Mdm) into SFR (M*) using 
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FIG. 1. — The distribution (shaded contours) of z ~ 2.2 halos in the dark 
matter accretion rate (left y-axis) versus halo mass plane. Also major merger 
fractions are displayed (red contours). Associated SFRs are indicated (right 
y-axis) assuming an effective star formation efficiency e = 1 in eq. ff] SINS 
galaxies at 2 < z < 2.5 are indic ated as disks (circles) or mergers (stars) 
based on the Shapiro et al. (2008) classification, and others (points) not in- 
vestigated by Shapiro et al. (2008). Their SFRs are based on their Ha fluxes 
corrected for extinction using A[j a = 0.8. Their halo masses were determined 
by assuming that the observed d isk maximum rotation vel ocity is equal to the 
circular velocity of the halo (Forster Schreiber et al. 2006), and only galaxies 
with M hillo > 10 1L25 Mq are included. The SFRs for the SMGs are shown 
as the horizontal lines in the upper part of the figure (because the SMGs 
are compact, their halo masses cannot be reliably inferred from the observed 
gas motions). The left y-axis indicates what dark matter accretion rates are 
needed to account for the observed SFRs when assuming e = 1 . The arrow 
indicates the shift in the galaxy positions if e = 0.5. 



the baryonic fraction tjb = 0.18 and an effective star formation 
efficiency e, which is a free parameter used to interpret the 
results: 



M* = r/ B x exM c 



(2) 



In the 'cold flow' regime ( M ha io 5, lO l2 M r ,; 
Birnb oim & Dekell 120031; iKeres et all 120051; lOcvirk et all 
2008) eq. |2]is a plausible measure of the baryonic accretion 
rate. At larger masses the accretion rate is lower and is 
controlled by the cooling time in the hot virialised baryonic 
halo gas. At much smaller masses it is strongly reduced 
by outflows generated by supernovae feedback. In the cold 
flow regime, the cold gas (which may be clumpy) is fed at 
approximately virial velocity via filaments directly into the 
halo center, where it accumulates onto the galaxy. Moreover, 
considering the case where gas is stripped off incoming 
galaxies, our estimated merger duration is related to the 
dynamical time of the halo, on which the gas will fall to the 
central galaxy when the cooling time is short. For major 
mergers e may even exceed 1 , because the star formation burst 
they trigger can be shorter than the dark matt e r halo merger 
time scale (e.g. ISpringel & Hernquisl l2005t iTacconi et alj 
[20081) . 

3.1. SINS galaxies 

On Fig. Q] we overplot the SINS galaxies. Their SFRs 
are based on their Ha fluxes corrected for extinction using 



SFR [NT n yr ] assuming e±1 




540 



-500 



500 1000 1500 2000 
DM accretion rate [M yr -1 ] 



2500 3000 



FIG. 2. — The distribution of dark matter accretion rates for halos with 
masses (1.2 ±0.3) X 10 12 Mq alz~2.4 (corresponding to the halo mass and 
redshift of BzK- 15504), which are not undergoing a major merger. The SFRs 
on the upper axes are given assuming effective star formation efficiencies 
e = 1 and e = 0.5. The measured SFR of 140(-80,+110) M Q yr 1 in BzK- 
15504 is indicated by the asterisks and error bars. 



A-Ha = 0.8, and their halo masses were estimated by assum- 
ing that the observed disk maximum r otation velocity is equal 
to the circular velocity of the halo dForster Schreiber et alj 
2006). Fig. Q] shows that if e > 0.5 is assumed, the host ha- 
los of SINS galaxies with M ha]o > 10" 25 M lie in the region 
where most halos of their mass are expected to be concen- 
trated. Furthermore, for e > 0.5 the expected mass accretion 
rates are sufficient to account for the observed SFRs. Also, the 
predicted major merger fraction is small (< 0.5), consistent 
with observations. Although the statistic s are still small, we 
notice that the confirmed SINS mergers (ShapiroeLaD]2008; 
stars) have higher specific dark matter accretion rates than the 
confirmed disks {open circles), and therefore come from a re- 
gion where the halo major merger fraction is higher. 

The computed number density of halos with M « 10 1L5 - 
10 12 M Q is a few times higher than the observed number 
density of the galaxies the SINS sample is drawn from (Ta- 
ble [T). Possibly, the observed galaxies have typically high 
Mg a i/Mh a i , with other halos of comparable mass hosting 
fainter undetected galaxies. Also, some of the halos in this 
mass range may have already developed virial shocks that 
quench star formation. 

In should be noted that our estimated major merger dura- 
tion Tinker equals w 350- 1000 Myr at z w 2.2, depending on 
the mass ratio. This is very similar to the "observa ble" galaxy 
merger timescale often found in the literature (e.g. Conselice 
120061; lLotz et al. 2008). Therefore, the approximation taken 
here, i.e. considering directly the halo merger fraction and in- 
ferring the galaxy merger fraction from it, is probably a rea- 
sonable one with respect to the deduced merger fraction of the 
S INS galaxies. 

iGenzel et all d2006l) studied BzK- 15504 with high resolu- 
tion using adaptive optics, and concluded that it was a large 
proto-disk with no sign of a recent/on going major merger, and 
a SFR of 140(-80 ,+l 10) M^ vr" 1 . iForster Schreiber et alj 
(120061) . lLawet alJ (120071) . iGenzelet alJ (12008b . G. Cresci 
et al. 



(2008, 



in preparation ), 
Ivan Starkenburg et alj J2008) 



Bournaud et alJ d2008) and 
have found similar systems. 



Fig.|2]shows that for halos not undergoing major mergers, and 
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with masses equal to the halo mass of BzK- 15504, the typical 
dark matter accretion rate is 450 Mq yr" 1 , i.e. the typical 
SFR assuming e= 1 is « 80 M yr" 1 . About 15% of such ha- 
los have SFRs exceeding the 140 M Q yr" 1 observed in BzK- 
15504 (again assuming e = 1). Thus, the implied dark matter 
accretion rate in BzK- 15504 may be quite typical. Consid- 
ering the uncertainty in the measured SFR, the implied dark 
matter accretion rate is consistent with theoretical expecta- 
tions for e as low as w 0.5. 

We conclude that high star formation rates and large abun- 
dances of non-major merger, massive disks at z ~ 2 are con- 
sistent with expectations from ACDM simulations if accretion 
is in the 'cold flow' regime and the star formation efficiency 
is high. 

3.2. SMGs 

For ewl, the observed SFRs of the SMGs imply dark mat- 
ter accretion rates of « 2500- 6000 M yr" 1 (Fig. IB- When 
examining halos with accretion rates in this range that are un- 
dergoing major mergers, we find that their masses lie mostly 
in the range (2-6) x 1O 12 M and obey a log-normal distri- 
bution with a mean w 3 x 1O 12 M and a sa 0.25dex. We 
also find that their number density is w 5 x 10" 5 Mpc" 3 . This 
is only slightly larger than the observe d SMG densi t y (Ta- 
bleQ]), and supports the conclusions of iTacconi et al.l (2006, 
2008) that the SMGs represent major mergers. If the SMG 
phase is shorter than the halo merger duration, such that e > 1 , 
the implied number density is not much altered, but lower 
halo masses are found. E.g., if the SMG phase lasts only 
lOOMyr, the mean halo mass is w 1O 12 M , in which case 
SMGs could be members of the UWoptically selected galaxy 
populations that have recently experienced a dissipative major 
merger. The observed rotation velocities of SMGs are larger 
than the expected circular velocities of halos with these in- 
ferred masses. This is consistent with the SMGs being con- 
centrated major mergers where the rotation velocities peak 
close to the center. 

4. FATE AT z = 

4. 1 . Fate of SINS galaxies 

Fig. [3] summarizes the major merger history of halos from 
z ~ 2.2 to z = 0. For halos with initial masses typical of the 
SINS galaxies' halos, sa 40% will be accreted via minor merg- 
ers by more massive halos (representing groups or clusters) 
with final masses 10 13 M < M < 10 15 M at z = 0. Around 
one half of those halos will merge fully with the central sub- 
halo of the group/cluster, and the other half will remain satel- 
lite subhalos. The other « 60% remain "main branch" halos 
to z = 0. Of these, « 2/3 undergo at least one major merger 
during their evolution to z = 0. Their final halo masses are 
10 12 - 3 M < M < 10 13 3 M , so their associate d galax ies may 
become massive ellipticals (cf. IConroy et al.l ((2008)). The 
other w 1/3 do not undergo any future major mergers, and 
grow to a mass 10 1L8 M < M < 10 12 5 M at z = 0. Thus, 
these may evolve via secular evolution into bulges and later 
possibly grow a new disk. 

4.2. Fate of SMGs 

A popular scenario is that the large central mass densities 
of SMGs and of z ~ 2 compact red sequence gal axies are re- 
duced by z = via dry dissipationless merg ers dvan Dokkum 
120051; iBell et all l2006t iNaab et all 12006b . ITacconi et al] 
(2008) show (in their Fig. 5) that this requires that the SMG 
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FIG. 3. — The number of future major mergers that halos with initial 
(z PS 2.2) masses 10 1L5 M Q < M < 10 12 ' 5 M Q (top) and 10 12 3 M < M < 
10 12 8 Mq (bottom) undergo from z Ri 2.2 to z = 0. We distinguish between 
(a) "main branch" halos that undergo only minor mergers to z. = (filled 
blue column at 0), (b) "main branch" halos that undergo at least one ma- 
jor merger (filled blue columns at > 0) and (c) "disrupted" halos (red). These 
are halos that at some stage are accreted onto a more massive halo in a mi- 
nor merger event. For category (c), major mergers are counted only prior 
to merging with the larger halo. Around 60% (ss 70%) of the halos in the 
lower (higher) mass bin remain "main branch" halos. For both mass bins 
the mean number of future major mergers is f» 1, while fewer than 10% un- 
dergo more than 2 major mergers. These results do not depend on whether 
the initial halo at z ~ 2.2 is identified as a major merger or on its dark mat- 
ter accretion rate. We find that the mean number of major mergers a halo 
with mass M undergoes between z; and z/ < Zi is well approximated by: 



AW,(Zi,Z/,M)«0.13x(log( 



I0'°M 



)+l)(z,-z/). 



masses grow by about an order of magnitude by z = 0, assum- 
ing dry mergers w ith structurally similar systems, following 
iNipoti et all d2003l) . We find that of the halos we have iden- 
tified with the SMGs in jj3.21 w 70% remain "main branch" 
halos to z = (Fig. [3]), and that their masses grow by factors 
of m 3 — 30. This mass growth appears consistent with the re- 
quirement of the dry merger hypothesis. However, most of 
the mass growth does not occur via major mergers, since typi- 
cally only w 1 major merger occurs per halo to z = 0, as shown 
by Fig. [3] 

Moreover, we show in Fig. H] that most of the mass 
growth is achieved via mergers less intense than 10 : 1, 
which is qualitatively consistent with the idea that the 
growth of massive galaxies is not dominated by major 
mergers (e .g. Haus man & OstrikerJfT978t iMaller et all 120061: 
Masjedi et al. 2008). This is especially true for the halos that 
do not undergo further major mergers until z = 0. Such halos 
tend to grow in mass only by a factor of sa 3, and gain > 0.7 of 
their new mass via mergers less intense than 10:1. Also, the 
galaxies themselves probably grow even less than their dark 
matter halos. It seems unlikely that the mass accreted via such 
small h alos can be sufficien tly gas poor for minor dry merg- 
ers (e.g. lBurkertet al.ll2007l) to be an important growth mech- 
anism (unless the galaxies are effectively stripped of their gas 
befor e merging with the descendant of the SMG (Naab et al. 
l2007h . 

A simpler and more likely explanation (ITacconi et al.l2008h 
is that the high SMG densities trace only the central starburst 
region, and exclude more extended fainter envelopes of pre- 
merger stars. At lower redshifts these stars would become 
visible and be of greater relative importance as the starburst 
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one major merger to z = (filled red). The results do not depend on whether 
the halo at z ss 2.2 is undergoing a major merger or on its accretion rate. 



fades, giving rise to larger h alf-light radii and small er inferred 
densities (cf. simulations by Hopkins et al. (2008b)). 

5. SUMMARY 

We have constructed new halo merger trees from the 
ACDM Millennium Simulation. Our trees account for merger 



durations, and we use them to identify halos that are under- 
going mergers and to extract dark matter accretion rates. We 
show that the high star formation rates observed in rotating 
disks at z ~ 2 are plausibly consistent with the dark matter ac- 
cretion rates expected for halos not undergoing major merg- 
ers. Given the measured star formation rates of SMGs, and 
the observationally supported assumption that they are under- 
going major mergers, we infer their likely halo masses. Major 
mergers can neither lead to the complete transformation of the 
z ~ 2 disks to z = ellipticals, nor to the disappearance of the 
very high density SMGs from z ~ 2 to z = 0. Therefore, sec- 
ular/internal processes are likely important in the evolution of 
these high-redshift populations to present time. 
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